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UNIVERSITY  OF  SOUTHERN  CALIFORNIA 

JOINT  SERVICES  ELECTRONICS  PROGRAM 
CONTRACT  NO.  F49620-85-C-0071 

RESEARCH  IN  ELECTRONICS 

DIRECTOR'S  OVERVIEW 

This  report  summarizes  the  progress  made  under  the  Joint  Services  Electronics 
Program  for  the  period  1  January  1986  through  31  December  1986.  It  is  the  second 
annual  progress  report  on  the  three-year  contract,  F49620-85-C-0071.  The  report 
covers  the  thirteen  research  projects  being  supported:  four  in  solid  state  electronics, 
four  in  quantum  electronics,  and  five  in  information  electronics. 

The  report  also  summarizes  the  significant  accomplishments  that  have  been 
achieved  during  this  research  period. 

SIGNIFICANT  ACCOMPLISHMENTS 

Fundamental  Limitations  on  the  Photorefractive  Recording  Process  -  Armand  R 
Tanguay,  Jr.  -  Research  Unit  SS2-3 

The  modeling  of  the  photoconductive  process,  when  combined  with  earlier  work 
on  diffraction  from  volume  phase  gratings,  has  extended  the  understanding  of  the 
photorefractive  recording  process  and  its  fundamental  limitations.  This  work  has  been 
applied  to  the  Photorefractive  Incoherent-to-Coherent  Optical  Connector. 

New  Organic  Materials  for  Rectifying  Organic  -  Semiconductor  Interfaces  -  S  R  Forrest 

-  Research  Unit  SS2-4 

A  broad  new  class  of  phthalocyanine  (Pc)  based  organic  materials  has  been  found 
to  form  high  rectifying  energy  barriers  when  deposited  on  to  inorganic  semiconductors. 
Successful  01  contacts  have  been  formed  between  p-Si  and  copper  Pc. 

High  Tq  Low  Threshold  Lasers  -  E.  Garmire  -  Research  Unit  QE2-1. 

High  To,  low  threshold  InGaAsP/InP  lasers  at  3  pm  have  been  successfully  grown. 
These  lasers  were  grown  on  p-type  InP  substrates. 
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SOLID  STATE  ELECTRONICS 

HETEROJUNCTION  MATERIALS  AND  DEVICES  EMPLOYING 
ULTRATHIN  LAYERS  GROWN  BY  METALORGANIC  CHEMICAL 

VAPOR  DEPOSITION 
(MOCVD) 

Research  Unit:  SS2-1 
P.D.  Dapkus 

For  the  Period:  1  January  1986  through  31  December  1986 

Progress 

The  reconfiguration  of  our  MOCVD  reactor  to  grow  ultrathin  layers  is  now 
complete.  Layers  as  thin  as  20 A  are  routinely  grown  in  the  AlGaAs/GaAs  materials 
system.  The  doping  of  n-  and  p-type  layers  throughout  the  whole  system  is  also  complete. 
We  have  used  diethylzinc  as  the  p-type  dopant  and  have  been  able  to  grow  layers  that 
range  in  doping  from  5x10*6  cm‘3  to  greater  than  1x10*9  Cm*3.  A  major  difficulty 
with  this  dopant  is  the  limited  dynamic  range  available  at  any  growth  temperature.  The 
dynamic  range  is  limited  by  the  sticking  coefficient  of  Zn  in  the  material  and  the  flow 
controller  ranges.  At  present  only  an  order  of  magnitude  dynamic  range  is  controllable  at 
most  growth  temperatures.  A  modification  of  the  system  to  include  multiple  sources  or 
multiple  mass  flow  controllers  is  being  contemplated.  A  similar  though  less  severe 
problem  is  experienced  in  the  n-type  doping  of  these  materials.  We  are  currently  using  a 
commercially  supplied  Si2H6  doping  source.  Unfortunately,  the  temperature  behavior 
of  the  doping  suggests  that  the  source  is  SiH4.  Decomposition  of  Si2H6  in  the  source 
bottle  is  promoted  by  the  residual  water  vapor  in  the  source.  In  spite  of  these  difficulties 
doping  over  the  range  1x10* 6  to  1x10*9  cm_3  is  achievable  in  the  appropriate 
temperature  ranges.  With  these  sources  calibrated  we  are  able  to  concentrate  on  the 
technical  challenges  in  growing  two  of  the  more  interesting  electronic  devices  that  can  be 
fabricated  by  MOCVD  -  the  selectively  doped  heterojunction  FET  and  the  heterojunction 
bipolar  transistor  (HBT). 

We  have  begun  to  fabricate  selectively  doped  heterojunctions  and  to  characterize 
their  electrical  characteristics  with  capacitance  measurements.  Normal  and  inverted 
structures  are  being  fabricated  in  a  configuration  that  allows  the  determination  of  the 
interface  charge  at  the  heterojunction  and  the  conduction  and  valence  band  energy 
offsets.  The  carrier  concentration  profile  of  an  inverted  heterojunction  structure  is 
shown  in  Fig.l.  This  sample  shows  an  interface  charge  density,  a  =  3x10*0  cm'A  We  are 
studying  the  effect  of  growth  conditions  and  doping  source  purity  on  the  measured 
interface  charge  density.  When  sufficiently  low  interface  charge  densities  are  achieved 
we  will  concentrate  on  the  fabrication  of  transport  samples  to  optimize  the  low 
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temperature  mobilities  of  the  two  dimensional  electron  gas  and  on  the  fabrication  SDHT 
devices. 

The  acquisition  of  high  purity  Ga  source  materials  has  allowed  us  to  fabricate 
materials  with  background  carrier  concentrations  as  low  as  3  xlO1^  cm'3.  The  highest 
mobility  we  have  achieved  is  8  x  10^  cm^/V-sec.  We  will  attempt  to  use  triethylgallium 
and  triethylaluminum  rather  than  the  methyl  compounds  to  improve  the  C  incorporation 
and  improve  the  purity  of  the  AlGaAs  and  GaAs.  Recent  work  at  GTE  laboratories  has 
shown  that  mobilities  in  GaAs  as  high  as  2  X  10^  cm2/  V-sec  are  possible  with  ethyl 
compounds  utilizing  low  pressure  growth  techniques. 

Modeling  of  HBT  structures  including  structural  and  fabrication  enhancements 
has  been  competed.  We  have  developed  a  PC  based  Modeling  algorithm  that  employs  a 
lumped  circuit  model  of  the  transistor  in  which  the  circuit  element  values  are  based  on 
those  calculated  from  the  physical  structure  and  materials  properties  employed  in  the 
structure.  The  materials  structure  parameters  are  used  in  a  physical  simulator  to  estimate 
the  gain  in  the  structure.  An  enhancement  of  the  simulator  has  been  developed  that  allows 
us  to  calculate  the  transit  time  contributions  of  the  various  portions  of  the  transistor  these 
are  compare  to  the  values  calculated  from  traditional  formulae.  Fig  2  shows  the  output  of 
this  simulator  showing  the  distribution  of  the  transit  time  of  a  Pnp  transistor.  The  plot 
shows  both  the  electron  and  hole  transit  time.  Note  that  the  transit  ti,e  in  this  device  is 
dominated  by  the  time  to  cross  the  neutral  base  region.  SPICE  simulations  of  the 
transistor  have  also  been  carried  out  to  compare  with  the  lumped  model  and  a  model  that 
converges  to  the  results  of  SPICE  and  the  physical  simulator  in  the  appropriate  limits  has 
been  developed.  With  this  model  we  have  considered  the  performance  of  microwave  and 
switching  transistors  in  the  Npn  and  Pnp  configuration  with  various  fabrication  and 
structural  enhancements.  The  results  of  these  simulations  suggest  that  the  major  limiting 
factor  in  the  high  speed  performance  of  HBT's  is  the  RC  time  constants  associated  with 
the  various  layers  and  junctions.  Pnp  transistors  can  be  designed  that  approach  the  speed 
of  the  fastest  Npn  designs.  Little  improvement  in  the  speed  performance  can  be  achieved 
by  incorporating  graded  base  layers  since  the  performance  is  not  limited  by  the  base 
transit  time.  The  major  improvement  can  be  had  by  introducing  a  self-aligned  structure 
and  by  reducing  extrinsic  capacitances  in  the  device  design. 

We  have  begun  to  fabricate  Pnp  transistor  structures  to  verify  the  gain 
characteristics.  The  early  structures  will  not  be  constructed  with  geometries  that  are 
conducive  to  high  speed  operation,  but  rather  will  concentrate  on  the  achievement  of  high 
gain. 


D.  A.  Sunderland  and  P.  D.  Dapkus,  IEEE  Electron  Device  Lett.,  EDL6,  648 
(1985). 

S.  P.  Den  Baars,  B.  Y  Maa,  P.  D.  Dapkus,  A.  D.  Danner  and  H.  C.  Lee,  J.  Crystal 
Growth,  77,  188  (1986). 

D.  A.  Sunderland  and  P.  D.  Dapkus,  IEEE  Trans.  Electron  Dev.  (Feb.  1987). 

H.C.  Lee,  A.  Hariz,  P.  D.  Dapkus,  A.  Kost,  M.  Kawase,  and  E.  Garmire,  submitted 
for  publication  in  Appl.  Phys.  Lett. 

List  of  Presentations 

D.  A  Sunderland  and  P.  D.  Dapkus,  Talk  presented  at  1986  Device  Research 
Conference,  Amherst,  MA. 

S.  P.  Den  Baars,  B.  Y.  Maa,  P.  D.  Dapkus,  A.  D.  Danner,  and  H.  C.  Lee,  Talk 
presented  at  ICMOVPEIII ,  April  1986,  Los  Angeles,  CA. 

H.C.  Lee,  A.  Hariz,  P.  D.  Dapkus,  A.  Kost,  M.  Kawase,  and  E.  Garmire,  Talk 
accepted  for  presentation  at  the  Topical  Meeting  on  Photonic  Switching,  Lake 
Tahoe,  March,  1987. 
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SOME  INVESTIGATIONS  OF  THE 
KINETICS  AND  MECHANISM 
OF  MOLECULAR  BEAM  EPITAXIAL  GROWTH 

Research  Unit:  SS2-2 
A.  MADHUKAR 

For  the  Period:  1  January  1986  through  31  December  1986 
Research  Objective 

The  task  proposed  under  work  unit  SS4-2  defined  its  objective  to  be  (a)  systematic 
growth  of  GaAs/AlxGai.xAs(100)  quantum  well  structures  grown  under  condition 

suggested  by  our  earlier  reflection-high-energy-electron-diffraction  (RHEED)  intensity 
dynamics  studies  to  be  optimized  and  (b)  to  continue  exploring  influence  of  growth 
conditions  dependent  surface  kinetic  effects  via  examination  of  RHEED  intensity 
behavior.  The  following  describes  the  progress  towards  these  objectives  during  this 
reporting  period. 

Progress 

(a)  During  this  period,  a  major  step  forward  was  the  identification  by  us  of  the 
presence  of  irreversible  behavior  of  the  static  and  steady  state  RHEED  specular  beam 
intensity  and  its  implications  for  the  existence  of  metastable  surface  step  density  which  is 
smoother  than  the  stable  static  surface  behavior.  This  is  shown  in  fig.  1 . 

(b)  The  above  noted  observation  of  an  intrinsic  behavior  of  the  surfaces,  for  the 
first  time,  allows  reproducibility  of  growth  conditions  from  day  to  day  on  the  same  MBE 
machine  and  from  one  machine  to  another,  independent  of  the  uncertainties  in  the 
customary  reliance  on  thermocouples,  infrared  pyrometers  and  ion  gauges  for 
determination  of  substrate  temperature  (Ts)  and  group  V  pressure  (Pas4),  respectively. 

(c)  A  series  of  single  quantum  well  structures  with  well  width  varying  between  3 
monolayers  (ML)  and  40  ML  were  grown  under  RHEED  determined  identical  growth 
conditions.  Equally  significant,  the  growth  condition  chosen  corresponds  to  the  optimal 
conditions  -  the  metastable  surface  noted  under  (a)  above  for  growth  with  the  new 
approach  of  growth  interruption,  and  the  maxima  in  the  steady  state  intensity  in  the  (Ts, 
Pas4)  two  dimensional  growth  parameter  space  at  the  chosen  growth  rate  for  the 
customary  practice  of  growth  without  interruption. 

Under  separate  (AFOSR)  sponsorship  the  photoluminescence  (PL)  and  PL 
excitation  (PLE)  behavior  of  these  samples  was  systematically  examined.  The  observed 
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PL  linewidth  V  of  the  first  electron  to  heavy  hole  free-exciton  recombination,  as  a 
function  of  the  well  width  (dw),  for  wells  grown  without  interruption  is  shown  in  fig.  2. 
The  nearly  T  ~  dw-l  behavior  demonstrated  for  the  first  time  invalidates  the  customarily 
used  notion  and  model  of  fluctuations  in  the  well  width  being  the  dominant  scattering 
mechanism  in  high  quality  samples.  Rather,  it  is  the  fluctuation  in  the  interfacia] 
chemical  composition  giving  rise  to  fluctuations  in  the  band-edge  discontinuity  on  the 
scale  of  the  exciton  size  and  a  shorter-ranged  alloy  disorder  scattering  which  control  the 
PL  linewidth.  The  solid  curves  in  fig.  2  show  the  behavior  of  T  calculated  in  a  theory 

based  upon  this  view  of  the  nature  of  high  quality  interfaces.  We  parenthetically  note  that 
the  absolute  values  of  T  seen  in  fig.  2  are  amongst  the  narrowest  line  widths  ever  reported 

for  quantum  wells  in  this  thickness  regime  and  testify  to  the  value  of  RHEED  studies  in 
arriving  at  optimized  growth  conditions. 

List  of  Publications 

1.  P.  Chen,  J.Y.  Kim,  A.  Madhukar  and  N.M.  Cho,  Jour.  Vac.  Sc.  Tech.  B4,890 
(1986). 


List  of  Presentations 

2.  F.  Voillot,  J.Y.  Kim,  W.C.  Tang,  A.  Madhukar  and  P.  Chen,  paper  presented  at  2nd 
International  Conference  on  Superlattices,  Goteborg,  Sweden  (Aug.  11-15,  1986). 
Proceedings  in  Superlattices  and  Microstructures  (In  Press). 

3.  J.Y.  Kim,  F.  Voillot,  P.  Chen,  A.  Madhukar  and  N.M.  Cho,  paper  presented  at  27th 
Electronic  Materials  Conference,  Amherst,  Mass.  (June  24-26,  1986);  Jour.  Elec. 
Materials,  15,  317  (1986). 


PL  LINEWIDTH  (meV) 


ELECTROOPTIC  MATERIALS  FOR  OPTICAL  PROCESSING  AND 

COMPUTING  DEVICES 


Armand  R.  Tanguav,  Jr. 
Research  Unit  SS2-3 

Progress 


During  the  most  recent  contract  period,  significant  progress  has  been  achieved  in 
the  modeling  of  photoconductive  processes  in  media  such  as  bismuth  silicon  oxide  (BSO). 
Such  modeling,  when  combined  with  studies  of  the  polarization  properties  of  diffraction 
from  photorefractive  volume  phase  gratings  (as  were  developed  and  reported  during  the 
previous  contract  period)  significantly  extends  our  understanding  of  the  photorefractive 
recording  process  and  its  fundamental  limitations. 

Analytical  modeling  of  the  photoconductive  processes  has  emphasized  a 
perturbation  series  study  of  the  photorefractive  recording  of  multiple  image  structures. 
The  immediate  application  of  this  work  has  been  for  a  more  detailed  understanding  of  the 
Photorefractive  Incoherent-To-  Coherent  Optical  Converter  (PICOC),  a  novel  two 
dimensional  spatial  light  modulator  that  has  been  explored  under  this  contract.  However, 
the  resulting  perturbation  series  model  applies  to  a  much  broader  range  of  devices. 
Details  of  the  analysis  and  its  application  to  PICOC  performance  is  provided  in  JSEP  Pub. 
1. 


One  immediate  consequence  of  this  modeling  has  been  the  observation  of  an 
inherent  inefficiency  in  the  photorecording  process  when  a  high  bias  electric  field  is 
applied  to  heighten  the  photoconductive  sensitivity.  This  inefficiency  applies  to  broad 
classes  of  photorefractive  recording,  not  just  to  PICOC.  The  inefficiency  pertains  to  the 
placement  of  photoinduced  charges  to  create  a  prescribed  charge  profile.  One  interesting 
technique  to  minimize  this  inefficiency  is  to  use  one  form  of  the  PICOC  process  whereby 
a  single  spatial  frequency  grating  is  recorded  with  high  bias  fields,  and  then  selectively 
erased  by  the  image-bearing  light  beam  with  no  bias  field  applied.  This  concept  will  be 
further  explored  in  the  following  contract  period. 

The  analytical  model  of  the  photoconductive  mechanism  is  inherently  limited  in  its 
accuracy;  more  accurate  models  require  numerical  techniques.  Such  a  numerical  model 
has  been  developed  during  the  contract  period,  and  applied  to  a  study  of  the  relative 
merits  of  "running  grating"  and  "stationary  grating"  photorefractive  recording 
techniques.  The  "running  grating"  technique  appears  to  be  most  effective  for  low 
modulation  depths,  i.e.,  for  weak  information  signals. 
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The  previously  developed  models  of  the  optical  polarization  properties  of  light 
diffraction  from  photorefractive  gratings  (JSEP  Pubs.  2  and  3),  when  combined  with  the 
photoconductive  models,  enable  a  more  accurate  description  of  diffraction  efficiency  for 
running  gratings  in  BSO.  Because  of  the  presence  of  optical  activity,  the  coupling 
efficiency  degrades,  typically  by  as  much  as  a  factor  of  two.  This  partially  explains  why 
electrooptic  coefficients  observed  in  photorefractive  measurements  have  been 
significantly  lower  than  those  measured  by  other  techniques.  Details  have  been  given  in 
JSEPConf.  1. 

List  of  Publications: 

1.  J.  Yu,  D.  Psaltis,  A.  Marrakchi,  A.  R.  Tanguay,  Jr.,  and  R.  V.  Johnson,  "The 
Photorefractive  Incoherent-To-Coherent  Optical  Converter,"  in  Photorefractive 
Materials  and  Applications,  J.  P.  Huignard  and  P.  Gunter,  Eds.,  Springer-Verlag, 
New  York  (in  press). 

2.  A.  Marrakchi,  R.  V.  Johnson,  and  A.  R.  Tanguay,  Jr.,  "Polarization  Properties  of 
Photorefractive  Diffraction  in  Electrooptic  and  Optically  Active  Sillenite  Crystals 
(Bragg  Regime),"  J.  Opt.  Soc.  Am.  B2,  321  (1986). 

3.  R.  V.  Johnson  and  A.  R.  Tanguay,  Jr.,  "Optical  Beam  Propagation  Method  for 
Birefringent  Phase  Grating  Diffraction,"  Opt.  Eng.  25,  235-249  (1986).  (Invited 
manuscript  in  Special  Issue  on  Materials  and  Devices  for  Optical  Information 
Processing.) 


List  of  Presentations 

1.  A.  Marrakchi,  R.  V.  Johnson,  and  A.  R.  Tanguay,  Jr.,  "Polarization  Properties  of 
Two-Beam  Coupling  With  Running  Gratings  in  BSO,"  1986  Annual  Meeting  of  the 
Optical  Society  of  America,  Seattle,  WA,  Oct.  1986. 
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ELECTRONIC  SURFACE  ANALYSIS  USING  NON-DESTRUCTIVE 
ORGANIC-ON-INORGANIC  SEMICONDUCTOR  CONTACTS 

S.R.  Forrest 
Research  Unit  SS2-4 

For  the  period  1  January  1986  through  31  December  1986 

Progress 

The  work  included  in  this  new  unit  began  in  April,  1986.  During  the  time  period 
between  April  and  December,  several  important  new  results  were  obtained  which  should 
lead  to  a  deeper  understanding  of  organic -on-inorganic  (OI)  semiconductor  contact 
barriers.  During  the  First  period  of  this  work,  the  OI  deposition  process  consisting  of  a 
high  vacuum  (P  <  10~8  Torr  obtainable)  deposition  system  was  brought  into  operation. 
These  are  the  highest  vacuum  conditions  yet  achieved  in  the  fabrication  of  OI  diodes,  and 
will  be  instrumental  in  the  coming  year  in  determining  the  effects  of  exposure  to 
environmental  conditions  on  the  size  of  the  OI  interface  energy  barrier  and  density  of 
interface  states.  In  addition,  a  two-fold  organic  purification  system  was  constructed 
whereby  the  crystalline  organic  material  is  First  cleaned  by  a  6-day  solvent  extraction 
method,  followed  by  a  4  to  6  day  gradient  sublimation  process.  Via  these  means,  the 
highest  purity  materials  can  be  obtained  for  future  work  on  the  fundamental  nature  of  the 
OI  interface. 

The  most  interesting  discovery  of  the  last  six  months  is  that  a  broad  new  class  of 
organic  materials  have  been  found  to  form  high  rectifying  energy  barriers  when 
deposited  onto  inorganic  semiconductors.  In  past  work,  only  perylenetetracarboxylic 
dianhydride  (PTCDA),  naphthalenetetracarboxylic  dianhydride  (NTCDA)  and  related 
compounds  were  found  to  exhibit  the  effects  of  both  space  charge  injection  as  well  as 
formation  of  an  appreciable  energy  barrier  to  room  temperature  electron  transport  at  the 
OI  interface.  Heretofore,  attempts  at  observing  the  same  effects  in  the  phthalocyanine- 
based  materials  have  not  met  with  success.  Nevertheless,  the  phthalocyanines  (Pc)  have 
excellent  space-charge  transport  properties,  and  form  extensive  crystalline  stacks  which 
suggest  that  they  should  also  exhibit  OI-rectification  properties.  In  addition,  the  Pc's  can 
be  made  to  exhibit  a  broad  range  of  optical  and  electronic  properties  simply  by 
substituting  various  metal  ligand  atoms  in  the  center  of  the  molecule.  For  these  reasons, 
branching  out  into  the  Pc's  is  of  central  importance  to  both  our  understanding  and  our  use 
of  the  OI  interface.  This  year  we  were  successful  in  forming  OI  contacts  between  p-Si  and 
copper  Pc  (CuPc).  The  barrier  height  measured  was  0.67  eV  (c.f.  0.75  eV  for 
PTCDA/p-Si  contacts).  The  variation  of  barrier  height  with  organic  composition 
suggests  that  the  mechanism  for  barrier  formation  is  an  intrinsic  property  of  the  organic 
thin  film.  An  additional  observation  is  that  the  diode  n-value  extracted  from  the 
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forward-biased  I-V  characteristics  is  nearly  double  that  of  the  PTCDA/p-Si  value  of 
n=1.5.  From  earlier  analysis  of  the  interpretation  of  the  01  n-value,  the  discrepancy 
indicates  the  existence  of  a  considerably  increased  surface  state  density  for  the  CuPc 
contact,  which  is  probably  also  responsible  for  the  lower  energy  barrier.  This 
observation  is  corroborated  by  distortion  observed  in  the  C-V  data  obtained  for  the  CuPc 
diode  as  opposed  to  the  PTCDA  device  (see  the  figure  below). 

Systematic  analysis  of  these  observations  are  now  underway.  Indeed,  we  expect 
that  a  detailed  experimental  probing  of  the  comparison  of  01  interfacial  properties 
observed  for  PTCDA,  NTCDA  and  CuPc  will  greatly  clarify  the  important  factors 
involved  in  01  energy  barrier  formation.  There  is  a  considerable  body  of  knowledge 
which  has  thus  far  been  accumulated  on  the  crystalline  structure  and  transport 
mechanisms  of  these  three  materials.  These  protypes  therefore  form  the  basis  of  our 
approach  to  understanding  01  heterojunctions. 

Progress  has  also  been  made  on  several  other  aspects  of  this  project.  For  example, 
in  past  work  it  has  been  shown  that  analysis  of  the  quasi-static  C-V  and  G-V 
characteristics  of  the  01  diode  gives  quantitative  information  about  the  density  and 
distribution  of  defect  states  residing  at  the  01  interface.  This  so-called  semiconducting 
organic-on-inorganic  surface  analysis  spectroscopic  (SOISAS)  technique,  however,  has 
several  shortcomings  including  dependence  to  bias  sweep  rate  and  low  sensitivity.  We 
have  recently  made  progress  in  developing  truly  a.c.  means  of  implementing  SOISAS. 
Initial  computer  simulations  indicate  a  greatly  enhanced  sensitivity  to  the  presence  of 
deep  levels.  Still  to  be  solved  is  a  means  of  simplifying  the  analysis  procedure  such  that 
the  data  can  be  interpreted  without  application  of  unnecessarily  complex  computer  data 
reduction  algorithms.  Finally,  in  anticipation  of  future  work  in  this  exciting  field,  we 
have  designed  an  apparatus  to  be  installed  in  our  high  vacuum  system  which  allows  for 
the  complete  fabrication  of  the  two  level  01  diode,  as  well  as  electronic  and  optical 
probing  without  breaking  vacuum  at  any  point  in  the  measurement  or  fabrication 
sequence.  This  capability  is  essential  if  we  are  to  determine  the  effects  of  environmental 
factors  (such  as  oxygen  and  water  adsorption  onto  the  inorganic  semiconductor  surface) 
on  the  formation  of  the  01  energy  barrier.  Additionally,  this  apparatus  can  be  used  to 
determine  the  utility  of  the  01  barrier  as  a  contaminant  gas  detector,  etc.  The  equipment 
necessary  to  implement  these  experiments  is  currently  on  order,  and  we  expect  to  begin 
assembly  by  early  January,  1987. 

List  of  Publications 

1.  S.R.  Forrest,  M.  L.  Kaplan,  and  P.  H.  Schmidt,  "Organic  thin  films  for 

semiconductor  wafer  diagnostics,"  invited  chapter  submitted  to  the  Annual  Reviews 
of  Materials  Science,  vol.  17,  1987,  Annual  Reviews,  Inc. 


List  of  Presentations 


1 .  S.R.  Forrest,  "Semiconductor  contacts  using  organic  thin  films,"  invited  paper  at 
Physics  and  Chemistry  of  Semiconductor  Interfaces  Meeting,  Salt  Lake  City, 
January  27-29,  1987. 
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QUANTUM  ELECTRONICS 

TOWARD  ROOM  TEMPERATURE  LASERS  IN  THE 
3pm  WAVELENGTH  REGION 

E.  Garmire 
Research  Unit  QE2-I 


For  the  period  1  January  1986  through  31  December  1986 

Progress 


First,  high  T0,  low  threshold  lasers  predicted  and  developed  in  the  previous  JSEP 
phase  were  Finally  successfully  grown  and  characterized,  using  a  new  p-DCC  geometry. 
They  fully  lived  up  to  their  predictions.  We  published  an  Applied  Physics  Letters  article, 
Tom  Hasenberg  completed  his  thesis  and  final  manuscripts  were  prepared  describing  the 
results.  A  patent  application  was  Filed. 

Second,  new  students  are  taking  data  on  the  high  temperature  operation  of  the  p- 
DCC  lasers,  which  should  lead  to  another  paper  when  Finished. 

Third,  the  new  students  have  done  a  few  practice  LPE  runs.  However,  the 
equipment  went  down  and  replacement  parts  have  taken  a  long  time  in  arriving.  As  a 
result,  no  growths  of  3  pm  material  were  completed.  The  LPE  system  should  be 

operational  by  January. 


List  of  Publications 

1.  T.  C.  Hasenberg  and  E.  Garmire,  "Low  Threshold  High  T0  InGaAsP/InP  1.3pm 
Lasers  grown  on  p-type  InP  Substrates",  Appl.  Phys.  Lett.,  49,  400  (1986). 

2.  T.C.  Hasenberg  and  E.  Garmire  "An  Improved  Au  Be  Contact  to  p-type  InP". 
Journal  of  Applied  Physics.  September,  1986. 

3.  T.  C.  Hasenberg,"Low  threshold.  High  T0  InGaAsP/lnP  1.3  pm  Lasers  Grown  on 
p-type  InP  Substrates  with  a  three-Melt  Technique"  Ph.D.  Dissertation  in  Materials 
Science,  May  1986. 

4.  T.  C.  Hasenberg  and  E.  Garmire  Characteristics  of  InGaAsP/InP  p-DCC  lasers 
grown  by  the  Three-Melt  Technique  Submitted  to  Journal  of  Quantum  Electronics, 
November  1986 


5.  T.C.  Hasenberg  and  E.  Garmire,  "Heterostructure  Laser  and  Related  Method": 
Patent  Applied  for,  June,  1 986. 


List  of  Presentations 

1.  T.C.  Hasenberg  and  E.  Garmire,  "Low  Threshold,  High  To  GalnAsP/InP  Lasers"  , 
CLEO  '86,  post-deadline  paper. 

2.  T.C.  Hasenberg  and  Garmire,  "Low  Threshold,  High  To  1.3  pm  Lasers  Grown  by  a 
Three-Melt  Technique"  ,  Device  Research  Conference  June,  1986. 

3.  T.C.  Hasenberg  and  E.  Garmire,  "Performance  Characteristics  of  P-DCC  GalnAsP 
Lasers",  International  Semiconductor  Laser  Conference,  October  ,  1986. 


A  SPECTROSCOPIC  STUDY  OF  BASIC  PROCESSES  IN  ELECTRICALLY 

EXCITED  MATERIALS 

Martin  A.  Gundersen 
Research  Unit  QE2-2 

Report  Period:  1  January  1986  through  31  December  1986 

Progress 

New  spectroscopic  methods  have  been  developed  and  applied  to  the  study  of  high 
power  gas  phase  switches.  These  methods  now  provide  quantitative  information  about 
the  populations  of  excited  states  in  a  high  current  discharge,  data  for  time  dependent 
microscopic  models  of  high  power  switches,  and  other  information. 

Results  include: 

•  Time-resolved  absolute  measurements  of  Balmer  populations  in  atomic  hydrogen 
using  laser-induced  fluorescence. 

•  Application  of  laser-induced  fluorescence  to  the  measurement  of  Stark  broadening 
in  hydrogen,  thus  providing  time-resolved  electron  densities,  and  overcoming 
several  deficiencies  associated  with  other  methods. 


Development  of  a  time  dependent  model  of  the  recovery  phase  in  a  hydrogen 
thyratron. 
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•  Experimental  data  along  with  a  model  showing  three-body  recombination  to  be 
important  in  high  power  thyratron  recovery,  thus  pinpointing  the  role  of  this 
recovery,  thus  pinpointing  the  role  of  this  recovery  mechanism  for  the  first  time. 

The  spectroscopic  studies  effort  is  having  an  impact  on  the  development  of  high 
power  devices  in  the  following  way.  The  understanding  of  fundamental  physical 
processes  is  leading  to  an  understanding  of  high  power  device  limitation.  Further,  the 
models  and  data  can  be  applied  to  the  development  of  new  devices.  Thus,  it  should  be 
possible  to  develop  predictive  models  of  devices  that  are  based  on  new  ideas.  Examples 
would  include  most  devices  involving  electron  impact  ionization,  such  as  glow  discharge 
switches,  and  optogalvanic -effect  switches,  and  the  BLT. 

The  BLT  (back-light  thyratron)  is  a  development  that  is  related  to  and  is  partially 
an  outgrowth  of  this  research  effort.  It  is  a  new  switch  that  holds  promise  for  replacing 
thyratrons  in  some  applications.  The  BLT  has  operated  at  35  kV  standoff  voltage  and 
7000  A  peak  current.  It's  advantages  include  glow  discharge,  gridless  operation,  cold 
hollow  cathode,  high  peak  current,  high  current  rate  of  rise.  It  is  triggered  by  an 
unfocused  light  source. 


List  of  Publications 

1.  "Summary  of  the  workshop  for  research  issues  in  power  conditioning,"  M.A. 
Gundersen,  editor,  University  of  Southern  California  (1986). 

2.  "A  semi-empirical  formalism  for  the  calculation  of  deep  level  wavefunctions  in  k 
space,"  H-H.  Dai,  M.A.  Gundersen,  and  C.W.  Myles,  Phys.  Rev.  B.  33,  8234  (1986). 

3.  "Measurement  of  excited-state  densities  during  high-current  operation  of  a 
hydrogen  thyratron  using  laser-induced  fluorescence,"  D.A.  Erwin  and  M.A. 
Gundersen,  Appl.  Phys.  Lett.  48,  1773  (1986). 

4.  "Determination  of  electric  field  and  electron  temperature  in  the  positive  column  of  a 
high-power  hydrogen  thyratron  from  non-intrusive  measurements,"  D.A.  Erwin, 
J.A.  Kune  and  M.A.  Gundersen,  Appl.  Phys.  Lett.  48,  1727  (1986). 

5.  "A  low  pressure,  light  initiated,  glow  discharge  switch  for  high  power  applications," 
G.F.  Kirkman  and  M.A.  Gundersen,  Appl.  Phys.  Lett.  49,  494  (1986). 

6.  "Resonant  behavior  in  the  luminescence  processes  of  N-doped  GaP,"  H.H.  Dai,  M.A. 
Gundersen,  C.W.  Myles,  and  P.G.  Snyder,  in  preparation. 
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Vakuumelectronic,  NTG-Fachberichte  95,  94  (1986). 

8.  M.A.  Gundersen,  R.  DeWitt,  A.K.  Hyder,  C.R.  Jones,  J.A.  Kune,  M.J.  Kushner, 
E.P.  Muntz,  G.  Schaefer,  and  P.F.  Williams, "Research  issues  in  power 
conditioning,"  Proceedings  1986  Seventeenth  Power  Modulator  Symposium, 
Seattle,  WA,  June  23-25,  1986. 

9.  C.G.  Braun,  D.A.  Erwin,  G.F.  Kirkman,  and  M.A.  Gundersen,  "A  linear  thyratron 
for  developmental  research,"  Proceedings  1986  Seventeenth  Power  Modulator 
Symposium,  Seattle,  WA,  June  23-25,  1986. 


OPTICAL  SWITCHING 

William  H.  Steier 
Research  Unit  QE  2-3 

For  the  period  1  January  1986  through  31  December  1986 

Progress 

Photorefractive  materials  are  a  very  promising  nonlinear  medium  for  use  in 
optical  switching,  phase  conjugation  and  optical  information  processing.  The  high 
resistivity  III-V  semiconductors,  particularly  GaAs,  are  of  particular  interest  because  of 
their  fast  photorefractive  response  and  their  possible  use  in  the  infrared  at  diode  laser 
wavelengths. 

Two  wave  mixing  in  photorefractive  materials  can  be  used  for  optical  switching, 
optical  gain,  and  is  a  convenient  method  for  measuring  photorefractive  material 
parameters.  We  have  completed  a  series  of  experiments  on  GaAs:Cr  at  1.06  pwhich 
!  measured  the  two  wave  gain  as  a  function  of  grating  wavelength  and  applied  D.C. 

^  electric  Fields.  The  grating  formation  time  was  also  measured  as  a  function  of  applied 

1  Field  and  laser  intensity.  The  "moving  grating"  technique  was  also  investigated  to  further 

enhance  the  gain  when  applied  Fields  are  used. 

These  measurements  were  compared  to  theory  and  found  to  agree  reasonably  well 
|  when  the  correct  ratio  of  electronic  recombination  coefficient  to  electron  mobility  was 

'  assumed.  The  value  of  this  ratio  was  confirmed  in  work  supported  by  other  contracts 

'  which  measured  the  gain  using  r.f.  electric  Fields. 

t 


The  largest  gains  (~  1.0  cm-1)  were  observed  with  a  voltage  of  4.0  kV  across  the 
sample  and  the  grating  moving  at  optimum  velocity.  At  an  intensity  of  35  m  W/cm2,  a 
grating  formation  time  of  1.5  msec  was  measured.  The  formation  time  was  found  to  be 
inversely  proportional  to  the  intensity  so  that  a  formation  time  of  1 .0  psec  requires  ~50 
W/cm2.  With  applied  fields,  the  grating  wavelength  for  maximum  gain  is  2-3  (i  and  the 
gain  at  larger  wavelengths  (8-10  p.  is  an  order  of  larger  than  the  gain  with  no  applied 
field. 

Using  these  results,  the  potential  for  GaAs  in  several  different  optical  switching 
configurations  can  be  determined. 

List  of  Publications 

1.  G.  Albanese,  J.  Kumar,  and  W.  H.  Steier,  "Investigation  of  the  photorefractive 
behavior  of  chrome-doped  GaAs  by  using  two-beam  coupling,"  Optics  Lett.,  11, 
650  (1986). 

2.  J.  Kumar,  G.  Albanese,  W.H.  Steier,  and  M.  Ziari,  "Enhanced  two-beam  mixing 
gain  in  photorefractive  GaAs  using  alternating  electric  fields,"  to  be  published 
Optics  Lett.,  February,  1987. 

3.  J.  Kumar,  G.  Albanese,  and  W.H.  Steier,  "Photorefractive  two  -beam  coupling  with 
applied  radio  frequency  fields:  Theory  and  experiments.”  Submitted  to  JOSA  B. 

4.  J.  Kumar,  G.  Albanese,  and  W.H.  Steier,  "Measurement  of  two-wave  mixing  gain  in 
GaAs  with  an  applied  field  and  a  moving  grating"  in  preparation. 

List  of  Presentations 

1.  G.  Albanese,  J.  Kumar,  and  W.H.  Steier,  "Investigation  of  the  photorefractive 
behavior  of  chrome-doped  GaAs  using  two-beam  coupling,  "Optical  Society  of 
America,  Seattle,  WA,  October,  1986. 

2.  J.  Kumar.  G.  Albanese.  W.H.  Steier.  "Photorefractive  two-beam  coupling  in 
GaAs:Cr  with  applied  r.f.  fields  and  moving  grating,"  submitted  to  CLEO  '87, 
Baltimore,  MD. 
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SELF-PUMPED  OPTICAL  PHASE-CONJUGATING  LASER 
RESONATOR  CAVITIES 


Jack  Feinberg 
Research  Unit  QE2-4 

Report  Period:  1  January  1986  through  31  December  1986 

PROGRESS 

1)  We  demonstrated  that  two  separate  argon-ion  lasers  can  be  locked  together 
using  optical  phase  conjugation.  The  technique  uses  four-wave  mixing  in  a 
photorefractive  crystal  of  barium  titanate,  and  can  be  extended  to  more  than  two  lasers. 
The  barium  titanate  crystal  is  used  in  its  "self-pumped"  mode  to  optimize  the  gain  of  the 
system.  The  technique  takes  advantage  of  the  slow  response  time  t  (with  low-power 
illumination)  of  barium  titanate  to  keep  the  lasers  locked  to  within  a  frequency  of  1/ 1. 
The  two  argon-ion  lasers  were  phase-locked  to  much  less  than  1Hz  for  an  indefinite  time. 
The  extension  of  this  technique  to  locking  together  semiconductor  lasers  will  depend  on 
increasing  the  response  of  the  nonlinear  material  in  the  infrared  region  of  the  spectrum, 
and  can  be  accomplished  by  cooling  the  crystal  close  to  its  transition  temperature  at 
10°°C. 

The  results  of  this  work  were  published  in  Applied  Physics  Letters. 

2)  We  published  a  lengthy  study  on  the  origin  of  the  photorefractive  charge 
carriers  in  barium  titanate.  Our  experiments  show  that  oxygen  vacancies  play  an 
important  role  in  determining  the  effective  density  of  charge  carriers  in  barium  titanate. 
These  vacancies  may  form  complexes  with  iron  impurities  to  provide  the  donor  and 
acceptor  sites  required  for  photorefractive  activity.  Understanding  the  source  of  the 
charge  carriers  is  essential  for  further  experiments  to  alter  the  nonlinear  optical 
characteristics  of  this  material. 

This  paper  was  published  in  the  Journal  of  the  Optical  Society  of  America  B. 

3)  We  remeasured  the  Pockels  and  piezoelectric  coefficients  of  barium  titanate. 
Accurate  values  of  the  Pockels  coefficients  are  necessary  for  choosing  between 
competing  models  for  the  photorefractive  effect,  and,  in  particular,  for  determining  the 
magnitude  of  the  hole  and  electron  densities  in  these  crystals.  We  found  that  previous 
measurements  by  others  gave  significantly  smaller  values  for  the  Pockels  coefficients 
because  they  ignored  the  piezoelectric  effect  in  their  measurement.  Our  results  imply 
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that  there  is  more  competition  than  previously  suspected  between  holes  and  electrons  in 
barium  titanate. 

This  manuscript  is  in  preparation. 

4)  We  invented  a  new  device:  a  phase  conjugator  that  can  change  a  high-order, 
distorted  laser  mode  into  a  low-order,  undistorted  mode.  The  device  uses  a  "ring"  self- 
pumped  phase  conjugator  based  on  a  barium  titanate  crystal,  which  normally  generates 
the  phase-conjugate  replica  of  an  input  beam.  However,  we  found  that  by  placing  a 
Faraday  cell  in  the  ring  of  the  conjugator,  the  device  can  be  made  to  produce  a  different 
mode,  one  which  is  no  longer  the  phase-conjugate  mode.  The  Faraday  cell  breaks  the 
time-reversal  symmetry  of  the  ring,  and  makes  the  optical  path  length  different  for  light 
propagating  in  different  directions  around  the  ring.  The  barium  titanate  crystal 
compensates  for  for  the  asymmetry'  by  choosing  a  different  output  mode  that  has  a 
slightly  different  propagation  constant.  By  varying  the  asymmetry  caused  by  the  Faraday 
cell,  the  shape  of  the  output  mode  can  be  reproducibly  controlled. 

A  16mm  movie  showing  these  effects  was  made  by  the  principal  investigator.  This 
results  of  this  research  have  been  submitted  to  Optics  Letters. 
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INFORMATION  ELECTRONICS 

ANALYSIS  AND  SYNTHESIS  OF  PARALLEL  PROCESSING  SYSTEMS 

Dan  I.  Moldovan  &  G.  Bekey 
Research  Unit  IE2-1 

For  the  period  1  January  1986  through  31  December  1986 

Progress 

In  1986  the  research  has  focused  on  the  following  two  main  issues: 

Parallelism  analysis  in  rule-based  systems.  We  developed  techniques  for  reduction 
of  search  spaces  in  rule-based  systems  based  on  the  study  of  rule  inter-dependencies. 
Algorithms  for  allocation  of  rules  to  multiprocessors  were  studied.  Simulations  of 
parallel  executions  of  rule-based  systems  were  performed.  These  results  can  be  directly 
applied  for  speeding  up  the  execution  of  expert  systems.  The  following  papers  from  the 
1986  document  our  results:  Publications  2  and  5;  and  Presentations  2,  4,  5,  6,  and  7. 

Mapping  of  algorithms  into  array  processors.  This  is  a  continuation  of  our 
previous  work  in  the  area  of  mapping  algorithms  into  parallel  computers.  Specifically 
we  investigated  the  partitioning  of  algorithms,  algorithms  whose  data  communication  is 
described  by  permutations  and  tradeoffs  between  design  parameters  of  array  processors. 
These  results  were  presented  in  the  following  1986  publications  1,  3,  4,  and  6;  and 
presentations  1  and  3. 
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3.  D.  I.  Moldovan,  "ADVIS:  A  Software  Package  for  the  Design  of  Systolic  Arrays," 
accepted  to  IEEE  transactions  on  Computer-Aided  Design  of  Integrated  Circuits  and 
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Processors, "  submitted  to  IEEE  Transactions  on  Computers,  in  May  1985,  revised 
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5.  D.I.  Moldovan,  "Parallel  Processing  of  Production  Systems,"  submitted  to  IEEE 
trans.  on  systems,  Man  and  Cybernetics,  in  1986. 


6.  T.  C.  Lin  and  Dan  I.  Moldovan.  "Mapping  of  Algorithm  Permutation  into  Mesh- 
Connected  SIMD  Computers,:  submitted  to  IEEE  Transactions  on  Computers,  in 
1986. 
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BASIC  RESEARCH  IN  C^  DISTRIBUTED  DATABASES 


Victor  O.K.  Li 
Research  Unit  IE2-2 

For  the  period:  1  January  1986  through  31  December  1986 

Progress 

In  the  past  year,  we  have  focused  on  three  areas  of  basic  research  in  distributed 
databases,  namely,  performance  modeling  of  concurrency  control  algorithms,  query 
processing,  and  termination  protocols.  In  Pub.  1  and  Pres.  3,  we  have  developed  two 
performance  models  for  timestamp-ordering  algorithms.  We  believe  our  work 
represent  the  first  successful  attempt  to  model  timestamp-ordering  algorithms  without 
resorting  to  simulations.  The  model  for  locking  algorithms  that  we  have  developed  (Pub. 
5)  allows  one  to  include  the  effect  of  deadlocks.  The  Precedence-Assignment  model 
(Pres.  1)  is  a  unified  model  of  most  existing  concurrency  control  algorithms.  In  the  area 
of  query  processing,  we  have  developed  the  relation-partitioning  approach  to  query 
processing  (Pres.  4),  and  a  query  processing  algorithm  applicable  to  distributed  databases 
managed  on  a  local  area  network  (Pub.  3  and  4).  We  have  also  developed  termination 
protocols  in  distributed  databases  managed  on  networks  with  unreliable  components 
(Pres  2  and  Pub  6). 

Our  paper  describing  the  Event-Based  Reliability  Model  (EBRM)  for 
communications  networks  has  been  published.  This  paper  has  generated  much  interest  in 
the  research  community,  and  we  have  received  a  number  of  requests  for  reprints.  EBRM 
allows  one  to  model  and  evaluate  the  reliability/availability  of  a  network  in  which  the 
components  have  dependent  failures.  Details  of  this  model  are  given  in  Pub.  2. 
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IMAGE  TEXTURE  ANALYSIS  AND  RESTORATION 
USING  NONSTATIONARY  MODELS 

A. A.  Sawchuk 
Research  Unit  IE2-3 

For  the  period:  1  January  1986  through  31  December  1986 

Progress 

Work  on  multiple-resolution  texture  segmentation  has  continued  with  the  goal  of 
determining  the  best  techniques  for  utilizing  spatial  information  to  avoid  the  effects  of  the 
nonstationary  statistics  of  the  classification  features.  The  effort  concentrates  on  both  the 
low-level  (texture  feature  extraction  and  smoothing  of  feature  estimates)  and  at  the  high 
level  (segmentation  algorithms  incorporating  spatial  information). 

At  the  low  level,  have  concentrated  on  improved  methods  of  generating  useful 
spatial  information  to  use  in  the  final  classification  step.  We  are  examining  techniques  of 
low-level  texture  feature  smoothing,  such  as  the  use  of  adaptive  smoothing  filters 
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(previously  applied  only  to  gray  scale  intensity  values)  to  avoid  feature  mixing  at  region 
boundaries. 

At  higher  levels,  we  have  been  testing  a  variety  of  techniques  for  combining 
information  obtained  at  several  different  scales  of  spatial  resolution  in  order  to  find  the 
optimum  final  texture  classification. 

Several  different  methods  for  performing  the  final  hierarchical  classification  have 
been  tested.  The  methods  vary  in  complexity  and  effectiveness  and  the  current  goal  is  to 
select  the  method  that  performs  well  at  an  acceptable  computational  speed.  One  technique 
uses  the  coherence  or  compatibility  of  a  pixel  with  its  local  neighbors  to  refine  the 
segmentation  process,  especially  at  region  boundaries.  Iterative  relaxation  algorithms 
such  as  the*Rosenfeld-Hummel-Zucker  and  Peleg  methods  are  being  tried. 

Most  experimental  work  to  date  has  concentrated  on  the  segmentation  of  test 
images  consisting  of  a  mosaic  of  textures  arranged  as  rectangles  with  vertical  and 
horizontal  boundaries.  Our  current  work  is  aimed  at  examining  the  robustness  of  these 
methods  by  testing  them  on  more  general  types  of  texture  mosaic  patterns  having  rotated 
rectangular  regions  and  regions  with  non -rectangular  curved  boundaries. 
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REDUCED  MODELING  THROUGH  OPTIMAL  PHASE  MATCHING 

E.  A.  Jonckheere  and  L.  Silverman 
Research  Unit  IE2-4 

For  the  period:  1  January  1986  through  31  December  1986 


Progress 

LONG  TERM  OBJECTIVE  AND  MOTIVATION 

The  long  term  objective  of  this  research  is  to  derive  either  a  transfer  function 

model 

Pn*n  +  -  +  Po 

y(s)  =  - u(s) 

ansn  *  ...  +  oq 

or  a  state  space  model 
dx 

—  =  Ax  +  Bu,x  =  Cx  +  Du 

dt 

that  interpolates  an  experimental  frequency  response  data  {(u(/<w),y(/to))}. 

The  problem  formulation  grew  out  of  the  Investigator's  joint  work  with 
Livermore  Lab,  a  research  pertaining  to  the  control  of  SDI  laser  weapons.  The  finite 
elements  or  partial  differential  equation  modeling  of  very  large  space  flexible  structures 
is  most  likely  to  overwhelm  the  control  designer  with,  say,  10,000-20,000  vibration 
modes.  The  on-board  laser  beam  pointing  control  system  could  not  possibly  take  into 
consideration  all  structural  vibration  modes.  Besides  this  curse  of  dimensionality,  the 
higher  eigenfrequencies  could  be  as  much  as  300%-500%  in  error  and  the  modeling  of 
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structural  damping,  hysteresis,  etc.  will  be  extremely  difficult.  The  approach  taken  here 
radically  departs  from  modeling  from  "first  principles."  Here,  we  derive  a  model,  of 
acceptable  size  n,  that  matches  the  actuator-sensor  relationship  as  it  is  revealed  through 
on  line  experiments.  The  experimental  data  is  chosen  to  be  the  gain  and  phase  frequency 
response  data,  because  sine  wave  excitation  is  easy  to  reproduce  experimentally  and 
because  phase  and  gain  are  most  relevant  to  feedback  control. 

An  ideal  test-bed  appears  to  be  the  piezoelectric  beam  of  the  Aerospace 
Corporation  (Terry  Brennan).  It  is  a  flexible  beam  coated  with  a  piezoelectric  material 
that  acts  both  as  actuating  and  sensing  device.  The  rate  sensor  and  the  force  actuator  are 
co-located;  hence  the  transfer  function  is  stable/minimum  phase,  and  identification  is 
possible  from  the  phase  data  only.  Both  the  actuator  input  and  the  sensor  output  are 
electrical,  so  that  sine  wave  excitation  is  easy,  as  well  as  measurement  of  the  phase  shift 
between  u(ja)l  and  y(ju)l 

SPECIFIC  RESULTS  OF  THIS  YEAR 
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The  key  to  the  mathemancal  solution  of  the  above  formulated  problem  was  found 
in  the  course  of  this  JSEP  work  unit.  Consider  an  experimental  phase  data. 

\ij(ok) 

tfjCOlc)  =  ^  - 


uljcn) 


With  this  experimental  phase  data  we  construct  exp(2j  (wk)),  map  it  to  the  z- 
transform  domain,  and  an  FFT  yields  the  Laurent  expansion 


expl2j<p  (a)k))  =  ...  *  aj:  -  uq  +  a  .jz'1 
Then  we  construct  the  Hankel  matrix 


I  a. 2 

H=  |  a- 2  ... 
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The  reduced  order  model  (A.B.C.D)  that  matches  the  phase  data  is  derived  from 
the  singular  value  decomposition  of  H.  It  should  be  stressed  that  these  manipulations  are 
finite  dimensional  matrix  operations  and  no  complex  nonlinear  programming  is  required 
as  it  is  the  case  in  parametric  optimization. 


The  most  relevant  result  obtained  in  the  course  of  this  year  of  research  is  the 
following  bound  on  the  error  between  the  experimental  phase  data  <p  and  the  phase  <p  of  the 

model  of  order  n: 

v(0,/<p(j(0)  -  <p{j(0)l  < sin'U On+]  +On+2  +  ■■■) 

where  ai>C2>...>0  are  the  singular  values  of  the  Hankel  matrix.  The  salient  feature  is 
that  the  goodness  of  the  model  can  be  assessed  directly  from  the  experimental  data,  which 
is  not  the  case  for  optimal  parametric  modeling  through  nonlinear  programming. 
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RESEARCH  IN  COMPUTER  VISIONS 

Rama  Chellappa 
Research  Unit  IE2-5 

For  the  period  1  January  1986  through  31  December  1986 

Progress 

The  goal  of  this  research  unit  is  to  develop  model  based  algorithms  for  some 
computer  vision  problems.  We  have  made  significant  progress  over  the  last  twelve 
months  in  developing  algorithms  for  extracting  edges  from  noisy  images,  and  for 
extracting  linear  features  from  images  of  man  made  structures  such  as  airports  and 
buildings.  The  potential  application  of  this  research  is  in  smart  interpretation  of  aerial 
images. 

Both  these  algorithms  extend  the  edge  detector  developed  during  the  first  year  of 
the  project.  This  edge  detector  represents  the  pixels  in  the  noise  free  image  by  using  a 
space-variant  correlated  autoregressive  (AR)  random  field  model.  Due  to  this 
assumption,  the  first  and  second  directional  derivatives  are  approximated  as  a  function  of 


the  model  parameters  and  intensity  values  in  a  local  3x3  window.  An  edge  is  detected  if 
the  second  derivative  in  the  direction  of  the  estimated  maximum  gradient  is  negatively 
sloped  and  if  the  first  directional  derivative  and  a  local  estimate  of  variance  satisfy  some 
constraints. 

The  linear  feature  extractor  we  have  developed  consists  of  tracing  the  edge  points 
detected  as  mentioned  above,  linking  and  fitting  with  lines.  The  heuristics  used  in  tracing 
and  linking  are  specific  to  the  characteristics  of  the  the  AR  model  based  edge  detector. 
The  performance  of  the  linear  feature  extractor  is  illustrated  in  Figure  1 .  Edge  detection 
in  noisy  images  is  done  by  jointly  estimating  the  parameters  of  the  AR  model  and  the 
intensity  values  in  the  local  window  mentioned  above  using  a  reduced  update  Kalman 

filter  (RUKF).  The  overall  algorithm  is  interleaved  so  that  the  information  required  for 
edge  detection  at  (xq,\q)  is  computed  by  the  RUKF  centered  at  (aq  +  l,yo  +  1). 
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Research  L'nu  IE2-5  Figure  1  : 

Linear  segments  of  the  airport  image  extracted  using  our  method,  (a)  Original  airport 
image,  (b)  Outputs  of  the  AR  model  based  edge  detector,  (c)  Edge  diagram  indicating 
edge  pixels  as  dots  and  (d)  Linear  segments. 


